Understanding Jet Scaling and Jet Vetos in Higgs Searches 
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Jet counting and jet vetos are crucial analysis tools for many LHC searches. We can understand 
their properties from the distribution of the exclusive number of jets. LHC processes tend to show 
either a distinct staircase scaling or a Poisson scaling, depending on kinematic cuts. We illustrate 
our approach in a detailed study of jets in weak boson fusion Higgs production. 
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Higgs searches — at hadron colliders are typically 
plagued by large backgrounds. To extract Higgs signals 
we develop strategies to suppress backgrounds where the 
structure of QCD effects plays an important role. 

Searching for a Higgs boson produced we first reject 
backgrounds by reconstructing the invariant or trans- 
verse mass of the Higgs decay products. Based on the 
gluon-fusion production mechanism an additional dis- 
criminating feature turns out to be the boost of the Higgs 
boson, i.e. its recoil against a number of quarks or glu- 
ons. If the Higgs boson is usually produced with small 
transverse momentum we require a large opening angle 
of its decay products and veto hard jet activity. 

Combined Higgs search results based on exactly zero 
jets and one jet have recently been shown for the 
LHC [TJ |2] . For each case the kinematic extraction cuts 
can be optimized once we understand how signal and 
backgrounds fall into these jet bins. 

Extending these to two jets will include Higgs bosons 
produced in weak boson fusion (WBF), i.e. in associa- 
tion with two so-called tagging jets [3]. For this process 
the color structure of the signal together with the two 
forward jets forbids gluon exchange between the incom- 
ing quarks. Correspondingly, jet radiation into the cen- 
tral detector is suppressed. Requiring exactly two jets is 
equivalent to a central jet veto as a means to suppress 
backgrounds @]. 

As another example, the production of a Higgs boson 
in association with a W or Z boson [5] is mediated by an 
electroweak process. We therefore do not expect signif- 
icant jet activity in the signal, in contrast to Vy/Z+jets 
backgrounds. A jet veto makes use of this feature. 

In all these cases we need to compute and measure the 
numbers of jets in signal and background channels. From 
the theory perspective the requirement of observing ex- 
actly rijets jets is problematic [6 . Any computation based 
on parton densities which obey the DGLAP equation is 
per se jet inclusive, i.e. we always compute cross sec- 
tions including an unspecified number of collinear jets. 



In contrast, a jet veto corresponds to the probability 

l-Pvcto=^ = ^^, (1) 
Ctot Ctot 

in terms of the inclusive (tr n ) and exclusive (o~ n ) n-jet- 
associated cross sections and &o = <Ttot- 

The problem with Eq.Q is that we cannot really in- 
clude a sensitive detector region for the number of (visi- 
ble) jets. The collinear limit of the DGLAP equation is 
not a useful approximation because we only observe jets 
with finite transverse momentum. In an improved de- 
scription the initial state parton shower unfolds the par- 
ton splittings and produces more or less collinear initial 
state radiation. The approximate simulation of initial 
state radiation over the entire phase space critically lim- 
its our understanding of the exclusive nj e t s distribution 
and associated veto survival probabilities. 

Matrix element and parton shower matching is the key 
to simulating jet radiation at the LHC [7H3- It correctly 
describes the radiation of any number of jets over the 
entire radiation phase space, for example for pure QCD 
jet events |10| , W, Z or photon production in association 
with jets [TUlU3| . or Higgs production in gluon fusion [T5] 
as well as in weak boson fusion. 

We propose to study exclusive rij c ts distributions. A 
jet veto then becomes nothing but a cut on another well 
understood distribution. It will turn out that jet radia- 
tion is usually governed by one of two patterns, staircase 
scaling and Poisson scaling. We have access to both of 
them with current LHC sample sizes. Based on such scal- 
ing studies we can for the first time count exclusive jets 
in a completely testable framework. 

Staircase scaling — is defined as constant ratios of suc- 
cessive rtj ts rates for example for VF+jets or pure QCD 
jets production 

R(n+l)/n = " +1 = R ■ (2) 

The scaling parameter R depends on the core process 
and on the requirements on the jets, but not on n. This 
feature has been observed at UA2, Tevatron and the 
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LHC [mill]. Historically, Eq.Q is defined inclusively 
as R. However, exclusive and inclusive staircase scaling 
is equivalent with R = R [10] • Following the above ar- 
gument, we rely on the exclusive formulation to study 
the rijets distribution. From the constant jet ratio R we 
can derive the normalized distribution of the exclusive 
number of jets 
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At the LHC we can study staircase scaling for example 
in IF+jets and Z+jets production [TDHT21 [T5] . in pure 
QCD jet production [10], and in 7+jets production [15] . 
QCD production of Z+jets is defined by as Z radiation 
off strong jet production \M\ 2 cx aa". 

To observe staircase scaling it is crucial to consider 
total cross sections with as few kinematic cuts as possi- 
ble. While the reason for this scaling behavior from first 
principles is not entirely clear, we know that it is closely 
linked with the non-Abclian nature of QCD. In our sim- 
ulations we see that the majority of jets arise through 
initial state radiation (ISR) off the parton entering the 
hard process. This radiation mediates between the virtu- 
ality scales of a parton inside the proton and the parton 
entering the hard process. The large jet multiplicities 
constituting staircase scaling are then driven by further 
splitting of very few ISR quarks or gluons, i.e. we are 
sensitive to final state radiation patterns starting from 
hard ISR. This splitting ISR feature is driven by multi- 
ple gluon splitting, i.e. it only occurs for non-Abelian 
massless gauge bosons with a self coupling. 

In Fig. [l] we see that with an approximately constant R 
Sherpa [16] based on Ckkw matching [8] with at least 
four hard jets correctly reproduces experimental obser- 
vations. Note that the quoted values for (n + l)/n are 
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Figure 1: Simulated R( n +i)/ n distributions for Z+jets and 
i/+jets production via the effective g-g-H coupling. We only 
apply the basic cuts prj > 20(30) GeV and \yj\ < 4.5 for the 
solid (dotted) entries. The quoted n values are in addition to 
the two leading jets. 



counted in addition to the two tagging jets which define 
our core process. 

A proper analysis of the theoretical uncertainties shows 
that this scaling behavior is not affected jTIJt. Varying 
a a {mz) = 0.114-0.122 [T7] merely shifts R. A consistent 
variation of all renormalization, factorization and shower 
scales also leaves the scaling feature untouched, but with 
a large enough shift in R that such a shift should be 
viewed as a tuning parameter for jet merging |10] . Fixed- 
order QCD corrections leave the scaling untouched [18] , 
The only complication we observe in Fig. [I] is the varia- 
tion of the first entry i?i/o which corresponds to a need 
to properly define the hard process. Counting only jets 
above 30 GeV instead of 20 GeV alleviates this problem. 

In Fig. [2] we show the R( n +i)/n distribution for the 
WBF Higgs signal as well as the electroweak Z+jets 
channel at order \A4\ 2 oc a 3 a™. From signal and back- 
ground studies for WBF Higgs production [3J 0] we know 
that these processes have different jet radiation pat- 
terns. However, for the total cross section with mini- 
mal cuts we still observe an approximate staircase scaling 
R(n+V)/n — const. The slight drop in R( n +i)/n is due to 
different classes of Feynman diagrams contributing to the 
electroweak process at different jet multiplicities, includ- 
ing WBF topologies, Z bremsstrahlung, and WW/WZ 
pair production 

If we want to use jet scaling for Higgs searches, we need 
to consider background suppression cuts. As an illustra- 
tive example we study WBF Higgs production qq — > qqH 
with a decay H — ?> tt or H — > WW [3j. Fixed-order 
corrections to the inclusive rate are known to be moder- 
ate [TH] ■ We are only interested in the production process 
so neglect the Higgs decay products entirely. Starting 
from exactly two tagging jets defined as the hardest jets 
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Figure 2: -R( n +i)/n distributions for electroweak Z+]ets pro- 
duction and WBF Higgs production. Again, we only apply 
PTj > 20 GeV and \yj\ < 4.5 and count n in addition to the 
two leading jets. 
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Figure 3: R(„+i)/„ distributions for electroweak Z+jets pro- 
duction, and Higgs production in weak boson fusion. Unlike 
in Fig. [2] we now count only additional jets in the veto region 
denned by Eq.(|5|. 



with 

Ptj > 20 GeV 
2/12/2 < 



to I < 4.5 

2/2 1 > 4.4 
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> 600 GeV 



we compute the jet activity in the to-be-vetoed region 



p^ to > 20 GeV 



min 2/1,2 < 2/ voto < max 2/1,2 • (5) 



Instead of the two hardest jets fulfilling Eq.Q we can 
also use the most forward jets without our conclusions 
being affected. In Fig. [3] we first show the signal-like 
processes, now only counting jets inside the veto region 
Eq.([5]). These cuts reduce backgrounds while capturing 
the signal properties, most notably the large invariant 
mass of the tagging jets. Therefore, they do not have a 
major effect on weak boson fusion and electroweak Z pro- 
duction and the approximate staircase scaling persists. 
The suppression of the first additional jets merely reflects 
the suppression due to the color structure combined with 
the strictly enforced tagging jet structure. Slight devia- 
tions from a perfect staircase scaling are also expected 
because the WBF cuts in Eq.Q are given by the LHC 
analyses and not optimized to test scaling features. 

Poisson scaling — is a different scaling property which 
we observe for example in QCD Z+jets production after 
cuts. For a limited number of emissions up to n crlt it 
reads 
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in terms of n, the number of jets expected. Statistically, 
a Poisson distribution represents the number of positive 
outcomes of independent trials. We observe it in multiple 
soft photon radiation for example off an electron [2U] . Its 
theoretical derivation rests on two features: first, one 
splitting dominates, e.g. successive photon radiation off 
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Table I: Observables for staircase and Poisson scaling. The 
ratios R and R correspond to the exclusive and inclusive njets 
distribution, respectively. T(n, n) is the upper incomplete 
gamma function. In the last line we show the jet veto survival 
probability P vcto = a /a tot . 



an electron or gluon radiation off a quark; second, soft 
radiation is automatically ordered by the radiation angle, 
which means there is no combinatorial factor. 

From the previous discussion we know that such a radi- 
ation pattern does not occur for total jet rates in proton- 
proton collisions. The dominant splitting ISR via the 
gluon self-coupling instead leads to staircase scaling. 

This changes once we apply cuts forcing our events 
into a specific kinematic configuration. The form of the 
Poisson distribution is reminiscent of Sudakov factors, a 
solution to the DGLAP equation on which the partem 
shower approach is based. Such Sudakov factors model 
the leading non-splitting probabilities for example of par- 
tons entering the hard process coming from the proton. 
The approximation underlying the DGLAP equation is 
collinearity and a sizable difference in virtualities be- 
tween partons inside the proton and the hard process. It 
then corresponds to resumming the collinear logarithm 
arising from successive ISR. 

A similar kinematic situation we enforce through the 
WBF cuts Eq. Q : the partons have to generate the large 
invariant mass rrijj > 600 GeV which favors quarks in the 
initial state. These quarks can most efficiently evolve via 
successive soft-collinear gluon emission; i.e. they show 
the same radiation pattern as soft photons being radiated 
off a hard electron. The moment such a large (collinear) 
logarithm enhances jet radiation we observe Poisson scal- 
ing of the jet ratios in QCD Z+jets and -ff+jets, Fig. [4] 
Values of JJ^o > 1 reflect a much higher probability to 
observe jet radiation. 

Because Poisson scaling only enhances radiation as 
long as the large logarithm dominates the kinematics, 
for large rijets we do not approach zero, but the stair- 
case limit. Fitting Poisson curves to the first three bins 
in Fig. [i] gives us n = 1.4 for QCD Z production and 
n = 1.8 for Higgs production in gluon fusion. The last 
two bins follow a staircase pattern. A comparison with 
Fig. [2] confirms that we see consistent R values in the 
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Figure 4: R(„+\)/ n distributions for Z+jets production and 
Higgs production via the effective gluon-gluon-Higgs coupling. 
Unlike in Fig. [I] we now count only additional jets in the veto 
region defined by Eq.Q. The curves are fits to Eq.(|6|. 

staircase setup and the high-multiplicity Poisson regime. 

Comparing the QCD and electroweak Z+jets produc- 
tion processes we observe a very clear difference. While 
the staircase scaling property is indeed slightly sculpted 
due to the different subprocesses, the Poisson shape for 
the QCD process significantly enhances the first two jet 
radiations. 

Among other observables, Tab. |T] tells us how to com- 
pute a jet veto survival probability from the exclusive jet 
scaling. Based on our numerical results we cannot ex- 
pect to simply extract the two scaling parameters b and 
n and insert them into the resummed form for P V0 t o - O n 
the other hand, the distinct patterns shown in Tab. [T] al- 
low us to study the nj ets distributions before and after 
cuts, validate the appropriate simulation, and quantify 
the scaling ratios R( n +i)/n including their uncertainties. 
For WBF signal-like processes, where a measurement is 
difficult, it is crucial that we can understand and simulate 
the staircase scaling pattern which is typical for inclusive 
LHC processes. 

Outlook — In this Letter we have laid out a strategy to 
understand jet counting and the associated veto survival 
probabilities. The tool behind this study is the exclu- 
sive nj e t s distribution, as currently measured in several 
processes at the LHC. 

There exist two scaling types for exclusive jet rates 
at the LHC. First, staircase scaling for inclusive as well 
as for exclusive jet rates is defined by a constant ra- 
tio R = <j n+ i/a n . It occurs for total cross sections 
for example in W/Z/7+jets and pure QCD jet produc- 
tion [TnillSlEI] and is due to non-Abelian gluon splitting. 
In Higgs production staircase scaling is realized by all sig- 
nal and background inclusive rates and by the signal and 
electroweak Z production after tagging jet cuts. 

Second, Poisson scaling appears for processes in which 



jet radiation is enhanced by a large logarithm. It cor- 
responds to successive gauge boson radiation. In the 
past, this scaling has been used to describe all signal 
and background processes after tagging jet cuts [3J 2]. 
For QCD Z+jets and gluon-fusion Higgs production we 
confirm this result after requiring two tagging jets. 

The main effect of these two scaling patterns is visi- 
ble in the first emission. While for electroweak processes 
we find Ri/o = 0.15 — 0.3, typical QCD processes after 
WBF cuts range around R1/0 ~ 1-5 — 2.5. An exclu- 
sive jet requirement or central jet veto is a powerful tool 
to suppress backgrounds, provided that kinematic cuts 
drive the backgrounds into a Poisson regime and leave 
the signal of staircase type. 

An understood nj Cts distribution with and without 
Poisson-inducing cuts allows us to carefully study all as- 
pects of jet counting, including experimental and theoret- 
ical uncertainties. It should enable LHC searches which 
otherwise are plagued by severe theoretical uncertainties. 
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